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Development of an ellipsometry to the case of a coherent far infrared 
irradiation, low temperatures and small samples is described, including 
a decision of the direct and inverse problems of the convergent beam 
ellipsometry for an arbitrary wavelength, measurement technique and a 
compensating orientation of cryostat windows. Experimental results are 
presented: for a gold film and UBeis single crystal at room temperature 
(A=119 /um), temperature dependencies of the complex dielectric func- 
tion of SrTiOs (A=119, 84 and 28 //m) and of YBa2Cu307_5 ceramic 
(A=119 yum). 
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1. Introduction 

Ellipsometry is a well known technique^ of measurement of the op- 
tical constants of different substances. Despite of being known from the 
beginning of the century, ellipsometry has begun to develop intensively 
only together with lasers and computers. Work^ is one of the few early 
attempts to develop ellipsometry into far infrared. Recent years heavy 
Fermion systems and, especially, high-Tc superconductors has given a rise 
to the new efforts in the far infrared ellipsometry*^'^. In this paper we 
describe our method, ellipsometer and some experimental results. 
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2. Far-infrared ellipsometry on small samples 

2.1 Parallel or convergent beam? 

Far infrared is characterized by significant divergency of the beams due 
to the value of X/d ratio, where A is a wavelength and d is a beam diameter. 
It means, that it is impossible to create a thin parallel beam of far infrared 
as it is possible in visible. Thus, there are two alternative approaches to 
the far infrared ellipsometry on the small samples: 

1) the sample is in a wide parallel beam the electromagnetic field is of 
nonzero value on the edges of the sample, reflection and diffraction parts 
of the field are detected; 

2) the sample is placed in the focus of a lens or a mirror the field on the 
edges of the sample can be put equal to zero and diffraction effects may 
be neglected. 

Both these approaches suggest appropriate calculations of diffraction 
effects in the former and of the convergency of the beam in the latter case. 
In the case of a wide parallel beam diffraction field can be introduced in 
the parabolic-equation approximation of the scalar theory of diffraction 
as parabolic-like half-shadow zones from each edge of the sample^. At 
large angles of incidence (80°) these zones will be essentially inter-mixed 
just near the far edge of the sample, so the detected field will contain a 
diffraction part. Scalar theory of diffraction is not the sufficient one for 
the ellipsometry. This case the decision of the problem of a polarized 
plane wave diffraction on the impedance edge should be used. We suggest 
that such a decision can depend on the curvature radius of the edge as 
parameter, which is always unknown value in the experiment. That item 
needs a separate experimental study. Unfortunately, our sample holders 
ab initio were designed for the convergent beam geometry and do not allow 
us to carry out such a study. 

Is that diffraction effect from the edges a negligible one? To our opinion, 
that effect can be neglected only in the case of the low refiected samples. 
In the case of highly reflected samples the diffraction effect of the edges can 
be of the order of the small reflection effect to be measured. For example, 
gold mirror^'^ at A=100 /im and 80° due to reflection produces Aps=1.17° 
and ^'=44. 65°, where as usual p = tan ^' exp(iA). We have chosen a 
convergent beam (CB) ellipsometry. In the case of CB-ellipsometry we 
reduce to a minimum value the diffraction effects, but the problem of a 
correct account of the convergent beam reflection arises. We have found 
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Fig. 1. Optical scheme of the convergent-beam ellipsometer for far 

infrared: Li,L2,L^ — lenses with focal lengths /i,/2,/3, respectively; 
S — sample; Det — detector; tti — tts — fundamental planes; 9o — angle of 
incidence of the beam; the optimal position of polarizer and analyzer is 
between lenses Li and L2; and 1/a + = I//3. 



the decision of the direct problem of the CB ellipsometry, quite adequate 
to the experimental problems. 

2.1 Convergent beam ellipsometry 

The direct problem of the convergent beam ellipsometry is to calculate 
the detected signal in terms of the reflection matrix of the sample and the 
field distribution in the beam. Our decision^ is based on the lens's Fourier- 
transform capability^ . That decision simultaneously determines the optical 
scheme of the CB-ellipsometer (Fig. 1). Let's designate Exi{xi,yi) and 
Eyi{xi,yi) the complex amplitudes of x- and y-polarized beams in the 
input plane tti, then the field in the output plane 774 is: 



X4, 



J/4 



(x4,y4) = Rc(a;4,y4) 



Exi 
E, 



yi 



(a;4/i//2,-y4/i//2), 



where Rc is modified reflection matrix. It can be found in coordinates 
{9, (f), where 6 is the angle of incidence of a partial plane wave in its own 
plane of incidence and (p is the azimuth of this plane 



Rc{e,ip) = ALK{e,ip)AR{e,ip) 



(2) 



where 



R 



Rpp Rps 

Rsp Rss 



(3) 
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is a reflection matrix of a general kind^ for a partial plane wave, and 

1 



a 
-c 



-b 
d 



ad — be 



d 
-c 



-b 
a 



cos(0o) cos(0) cos(v9) + sin(6'o) sin(6') 
cos(0o) sin(99) 
cos(6') sm{{p) 
cos{ip) 



(4) 



(5) 



To calculate Rc - matrix in (x4, y4) coordinates, the variables are to be 
changed: 



cos{6) = < 3;4sin(^o) + cos(0o)\/ /I 



tan(99) = y4/ |sin(6'o) y^/l - xl-yj- cos(6'o)| 



(6) 



Placing an aperture in the output plane 114 we can filtrate some part of 
the output field, corresponding to the part of reflected plane waves. Lens 
L3 creates the image of the aperture in the detector plane tts. Detected 
signal in the ideal ellipsometer can be written as: 

Id{A) oc JJ^{Exi{x4,y4)cos{A) + Ey^{x4,y4)sm{A)}{. . .}* dx^dy^ (7) 

where A — angle of analyzer, {. . .}* means complex conjugation and the 
integral is to be taken over the aperture Ap in the plane ^4. To compare 
the detected signal in the case of the convergent beam to that one in the 
case of the alone plane wave expression should be rewritten: 



IniA) oc Ixx cos'^iA) + I^y sm{A) cos{A) + lyy sin^iA) 



(8) 



In the case of alone plane wave the detected signal has the same appearance 
as with the following differences: 



convergent beam: 



plane wave: 





= SApExiEl^dx4dy4 




Ixy 


= IlApiExiE*^ + E*^Ey^) dx4 dy4 


(9) 


^yy 


= IlApEy4Ey^dx4dy4 




Ixx 


= Ex^E*^ 




Ixy 


TP JP^ 1 77'* 77' 

- Ex^Ey^ + E^^Ey^ 


(10) 


^yy 


= Ey^E*^ 





4 



If a sample is an ideal mirror, expressions @ coincide with (10) ones. In 
the opposite case, each point of the output plane 7r4 has its own polarization 
ellipse and averaged coefficients (|9|) should be computed. 



In addition to the quantitative results above two qualitative notices 
can be done. 

a) The convergency of the beam will affect mainly on the averaged ellip- 
ticity of the beam. 

b) No matter in what point between the lenses Li and L2, together con- 
sisting a telescope, the reflecting sample is placed. 

2.2 Possibility of Multiple- Angle-of-Incidence (MAI) ellipsome- 
try 

Optical scheme of the CB-ellipsometer (Fig. 1) contains an interest- 
ing opportunity of MAI measurements. To carry out MAI measurements 
someone should move the output aperture Ap in the plane in x-direction 
perpendicularly to the optical axis. MAI ellipsometry can be used to de- 
termine the thickness of the film to be investigated. We checked such 
opportunity on pure substrate of SrTiOs. Some of results are in Fig. 2. 
Parameter x follows to the model computation line, while parameter 7 
deflects from it's line appreciably. Parameter 7 is mainly connected with 
a phase shift. It points out to the phase distribution in a real beam as a 
possible reason for that deflection. So, we can't use that simple scheme 
for MAI ellipsometry, although measurements at different central angles 
of incidence give more close results^. We suggest that measurements with 
uncoherent light can give better results. 

3. EUipsometric technique 

3.1 EUipsometric scheme 

We have chosen compensatorless scheme of the ellipsometer to simplify 
experiment, to work in a wide spectral region, and to avoid additional 
difficulties in solving the direct problem of CB-ellipsometry. We measure 
two parameters of the polarization ellipse - azimuth x and ellipticity 7 
- both averaged over the aperture. Experimentally measured values are 
analyzer angle Amin , at which the detected signal is of minimal value and 
relation of the minimal signal to the maximal one: 
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X = ^min (11) 
tan2(7) = lD{Amin)/lD{Amin^T^/'2) (12) 

These values can also be computed from (|^, (|9|). Equating the first deriva- 
tive of the detector signal (^) to zero, someone obtains: 

^exti = ^ arctan( 'j^ ) (13) 

Aext2 = Aextl ± 7r/2 (14) 

where A^xti and A(.xt2 are extrema of the function l£){A)^ one of those is 



minimum, another — maximum. In any case, to use (12) the function (^) 
should be computed at both points (|l3|) and (|T^). Computer program can 
readily distinguish which one is minimum and vise versa. Thus, we can as 
measure as calculate two parameters: averaged azimuth x and averaged 
ellipticity 7, that permits to decide the inverse problem. 



3.2 Inverse problem 

In brief, the numerical inversion problem is to find such values of pa- 
rameters of the reflecting system, that error function F has a minimum 
value^. For the multiple-angle-of-incidence (MAI) technique let us define 
the error function F as follows: 

M 

^ = E{(xr-xD' + (7r-7f)'} (15) 
1=1 

where M is the number of measurements (the number of incidence angles), 
xf^ and 7-" denotes angles from the ith measurement, Xi and 7^ are com- 
puted values. In this paper only a simplest reflecting system — an isotropic 
semi space — is discussed. In that case F is a function of two variables, 
e.g. n and k — refraction and extinction coefficients of the substance. It's 
necessary to study the surface of F{n, k) to draw the conclusion about 
existing of the single minimum of this function. 

Let's take SrTiOa single crystal as probe material, because of its large 
optical constants in far infrared, cubic structure and wide use as a substrate 
for YBaCuO films. Results of a simple computer experiment are in Fig. 3. 
Let's take the computed for the convergent beam values x^ and 7"* as 
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Fig. 3. Contour map of the 

error function i^(n, k) at n=22, 
K=13, A=119 /xm, M=l, 6*0=80°, 
aperture of diameter 30 mm, 
Li=L2=153 mm. 



Fig. 2. Ellipsometric parameters x 
and 7 vs a shift of a 6 x 40 mm sht 
in the plane tta of Fig. 1; crosses and 
circles — experiment, solid lines — 
computation, dashed line — intensity 
of the beam at s-polarization. 
measured ones, but x^l*^; i^) ^^id 7'^(n, k) we will calculate for the plane 

wave. This is a model of a convergent beam experiment with a plane wave 
treatment. Two conclusions can be drawn from Fig. 3: 

i) the error function F{n, k) has a single minimum, so the inverse prob- 
lem has a single decision; 

ii) the displacement of the minimum illustrates the influence of the 
focusing of the far infrared beam or difference between two approximations: 
of the plane wave and of the convergent beam. 

Comparing % = 11.18°, 7 = 7.02°, computed for the convergent beam, 
to X = 11-19°, 7 = 5.88°, computed for the plane wave, we can draw the 
third conclusion: 

iii) the convergency of a beam influences on value of 7 and practically 
does not influence on value of %• 



3.3 Problem of cryostat windows 

In a general case window is a thick anisotropic plate. In the case 
of a convergent beam windows must be thin. It means that 1-st, 2-nd, 
. . . , ith convergent beams, reflected inside window plate, all must have 
the same focal point. Because of this reason, due to its transparency in 
visible and far infrared and also due to its low temperature properties we 
use 20 jim mylar as material both for warm and cold windows. Mylar is 
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known to be an anisotropic material in far infrared . Rotating 20 /xm 
mylar foil between crossed polarizer and analyzer (PJ-A) we detected two 
approximately perpendicular axises of zero effect to the polarization of the 
beam. Maximum signal in that experiment was 8 x lO""^ of the signal 
in P\\A orientation at A=119 fim. Black polyethylene shows 1.4 x 10~^ 
at A=84 fim. Being mounted in an arbitrary orientation 4 windows and 
2 filters can give a significant polarization effect. We proposed following 
minimization and account of this effect. 

According to our development above, in presence of the windows, matrix 
^c{x4,yi) in equation (|^) should be modified: 

^c{x4,y4) ^T4{x4,y4)T3(x4,y4)Kc{x4,y4)T2{x4,y4)Ti(x4,y4) (16) 

where Ti(x4, 2/4) is 2 x 2 matrix of each window for every partial plane wave. 
But a strict account of the windows is practically impossible because we 
deal with the not normal incidence of a convergent beam onto a concave or 
convex anisotropic film with an unknown model of anisotropy. We prefer to 
minimize the windows anisotropy effect by their compensating orientation 
as it is described in Appendix and measure the rest of the effect due to the 
not normal incidence of the beam as a complex constant T\y, describing 
the transmission of all for windows. This value is used in solving of the 
inverse problem for the cold measurements. 



4. Experiment 

4.1 Setup 

According to the optical scheme of Fig. 1, far infrared ellipsometer was 
made up with water vapor electro-discharged laser being the source and 
liquid helium cooled photo-resistors Ge:Ga, Si:B, GaAS as the far infrared 
detectors. The section of the setup by plane of incidence is shown in Fig. 4c. 
All lenses are made up of polyethylene. Both polarizer and analyzer consist 
of two sheets: i) alumina film stripes on thephlon replica with .8 fim period 
and ii) metal mesh with periods 20x400 //m. Polarizers should be placed 
as close to the sample as possible to minimize polarization inhomogeneity 
of the beam due to that one of the optical elements, including polarizer 
and analyzer. The sample is in the atmosphere of a cryoagent. To filtrate 
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Fig. 4. Setup. 



the short wavelength irradiation a black polyethylene (in compensating 
orientation) surrounds the sample. 

Block-scheme of the setup is in Fig. 4a. By means of precision multi- 
revolution potentiometers P- and A-angles are transformed to the voltage. 
Accuracy and linearity of this transformation both are of order ±0.1°. To 
measure an intensity of light a lock-in technique at modulation frequency 
80 Hz is used. The polarizer is at fixed angle ±45°, the analyzer is rotated 
near the minimum of intensity and this part of the curve Id (^) is recorded 
in the computer file, and /©(^Imm) are calculated and analyzer 

begins to rotate to the A^ax angle. A linearity range of the detector does 
not exceed 10^, but we need to measure values of Imax/Imin of order 10^ 
for the metallic samples. So we use paper sheets as attenuators (each of 
order 3 at 119 /xm) to measure the signal at the same gain of the lock-in 
amplifier. Such measurement routine requires a time stability of the laser 
irradiation. 

The time instability of irradiation and the small movements of the 

sample in the cooling process (up to 0.3 mm) are two main sources of 
random errors. The whole cryostat is mounted on the movable optical 
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table to adjust the sample in x-direction. Errors due to tlie beginning of 
cooling process are good seen in the upper curves in Fig. 5. 

4.1.1 Room temperature results 

Some experimental results arc presented in the Table to illustrate the 
sensitivity of the ellipsomctcr. All these results were obtained at room 
temperature, at A=119 iJ,m without windows. Expected values for gold'''*^ 
are x = 45° — 0.36° and tan^(7) = 9.64 x 10"^. As can be seen from 
a transmittance measurements, accuracy of the ellipsometer is sufficient 
for gold's measurements. Increasing value of ellipticity 7, corresponding 
mainly to the phase shift, can be caused by uncertainties of the surface 
of the gold film. For UBeis we have computed Wp = 14600 cm~^ and 
Wt = 670 cm~^. Computed from these values free clcctrons's concentra- 
tion (m-e = 1) iVg = 2.35 X 10^^ is in good agreement with that one from 
Hall-effect measurements^^. 



Sample 


Measr. 




Polarizer 


angle, P 






values 


-45° 


0° 


+45° 


+90° 


No sample, 
transm. 


X 

tan^ (7) 


+0.13° 
5 X 10-6 


+0.25° 
5 X 10-6 


-0.10° 
5 X 10-6 


0° 

5 X 10-6 


Au 
film 


X 

tan^(7) 


+0.39° 
4.14 X 10"^ 


-0.01° 
9 X 10-6 


-0.60° 
6.28 X 10-^ 


0° 

7 X 10-6 


UBei3 
single cr. 


X 

tan^(7) 


+3.53° 
4.74 X 10-^ 




-3.41° 
4.53 X 10"^ 





4.2 Results of SrTiOs measurements 

Temperature dependencies of real and imaginary parts of the complex 
dielectric function of SrTiOs single crystal arc shown in Fig. 5. Such kind of 
behavior is caused by softening of low- frequency mode, centered between 80 
- 90 cm~^. We have measured different samples from independent sources 
and obtained different optical constants, especially at the frequency 84 
cm-^. Re-polishing of the samples does not change their optical constants. 
In literature^^-^^ low-frequency constants also different. We suggest that 
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Fig. 5. Temperature dependence of the 
complex dielectric function of SrTiOs 
single crystal at 3 wavelengths. 




Tc=91K 
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TEMPERATURE, K 

Fig. 6. Temperature dependence 
of the complex dielectric function 
of YBaCuO ceramic at wavelength 
119 /im. 



the lowest frequency mode is very sensitive to the quality of SrTiOa single 
crystals. 



4.3 Superconductive transition in YBaCuO ceramic 

Temperature dependence of the complex dielectric function of 1-2- 
3 ceramics at wavelength 119 /xm is shown in Fig. 6. Sintered on air 
and annealed in oxygen atmosphere sample was polished with diamond 
paste. The behavior of the complex dielectric function is determined by 
chaotically oriented small crystals with different optical constants for a—, 
6—, and c— axis. Porous structure of the surface of the ceramics also effects 
on the measured values of the complex dielectric function. 
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5. Appendix 
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Let 2x2 matrix T describes the transmittance of some window. Rota- 
tions of this window around the axes give the following results: 

^ (17) 
(18) 
(19) 

Transmission of the first two windows Wi and W2 can be written as: 
T12 = T2T1 



Til 


T12 




" hi 


h2 ' 


T2I 


T22 






t22 



Tuhl + Tl2i21 nih2 + Ti2i22 
T2lhl+T 22*21 T 21*12 + T 22*22 

If windows Wi and W2 are identical and W2 is rotated then 
T2 



■ 12 



Til 


T12 




*22 


±*21 


T2I 


T22 




_ ±*12 


*11 



(20) 



(21) 



and, putting the off-diagonal elements much less then the diagonal ones, 

*11*22 i *21 *22(*12 ± *2l) 
*ll(*21 =t *12) *11*22 =t *12 



*11*22 



1 (*12±*2l)All 
(*21 ± *12)A22 1 



(22) 



where '+' means the presence and ' — ' means the absence of {x —x) 
or (y — > —y) rotations. So the values of the off-diagonal elements can be 
reduced to the minimum. 

Let's now consider windows Wi and W2 as one window W12 and win- 
dows W3 and Wi as one window W34 (Fig. 3B). Then the reflection matrix 
R becomes: 



W34RW12 



Til 


T12 




T21 


T22 





Rpp Ftps 

^sp -^ss 





' hi 


*12 " 




*21 


*22 



RpphlT 11 Rpst22T 11 + Rpph2T 11 + Rsst22T 12 
RsphlT22 + -Rpp*llT21 + -Rss*2lT22 -Rss*22T22 



(23) 
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Due to the rotations (p^),([l^) matrix W34 should be changed according to 
and the final expression for the reflection matrix between two windows 
can be written as follows: 

W34R W12 = tut22 

Rpp {Rpst22 + Rpptl2 ^ Rsst2l) /hi ^24) 

(Rsptn ± Rpptl2 + Rsst2l)/t22 Rss 

where, as usual, '+' means the presence and '— ' means the absence of 
(x — > —x) or (y — > —y) rotations. So, practically full compensation of the 
cryostat windows effect can be reached. 
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